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Abstract. The 0-instanton solution of Painlcve I is a sequence (tt„^o) of complex numbers which ap- 
pears universally in many enumcrativo problems in algebraic geometry, graph theory, matrix models and 
2-dimensional quantum gravity. The asymptotics of the 0-instanton {u^fi) for large n were obtained by the 
third author using the Riemann-Hilbert approach. For fc = 0, 1, 2, . . . , the fc-instanton solution of Painleve 
I is a doubly-indexed sequence {u„ fc) of complex numbers that satisfies an explicit quadratic non-linear 
recursion relation. The goal of the paper is three-fold: (a) to compute the asymptotics of the 1-instanton 
sequence {u„^i) to all orders in 1/n by using the Riemann-Hilbert method, (b) to present formulas for the 
asymptotics of {u„ f^) for fixed k and to all orders in 1/n using resurgent analysis, and (c) to confirm numer- 
ically the predictions of resurgent analysis. We point out that the instanton solutions display a new type of 
Stokes behavior, induced from the tritronquee Painleve transcendents, and which we call the induced Stokes 
phenomenon. The asymptotics of the 2-instanton and beyond exhibits new phenomena not seen in and 
1-instantons, and their enumerative context is at present unknown. 
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1. Introduction 

1.1. The Painleve I equation and its 0-instanton solution. The Painleve I equation 

(1.1) -^u" + u^ = z 

is a non-linear differential equation with strong integrability properties that appears universally in various 
scaling problems; see e.g. [FIKN06]. The formal power series solution 

oo 

(1.2) uo{z) = z'^' ^n,oz-''''/' e z'/^C[[z-'/']] 

n=0 

and correspondingly the sequence (un,o) is the so-called 0-instanton solution of Painleve I. Substituting (1.2) 
into (1.1) collecting terms, and normalizing by setting uo,o = 1, implies that (unfi) satisfies the following 
quadratic recursion relation 

(l-o) Un,0 — Mri-1,0 ~ 2 ^ "<!,OMn-f,0, Uo,0 = 1- 

The 0-instanton solution (u„_o) of Painleve I is a sequence which plays a crucial role in many enumerative 
problems in algebraic geometry, graph theory, matrix models and 2-dimensional quantum gravity; see for 
example [DFGZJ95, GLM08, GMIO]. The leading asymptotics of the 0-instanton sequence (un.o) for large 
n was obtained by the third author using the Riemann-Hilbert approach; see [Kap04]. In [JKOl] (see also 
[GLM08, App.A]), asymptotics to all orders in n were obtained as follows 

(1.4) «„,o - A-2"+^r(2n - i) :^(l + V ^ j, n ^ 00. 

^ trnU(2n-i/2-fc)j 

In this expression, ui i are the coefficients of the 1-instanton series (defined below), A is the instanton action 

8V3 



(1.5) A 
and 

si 

(1.6) S, = -i-^. 
is a Stokes constant. 

1.2. The Painleve I equation and its fc-instanton solution. In this paper we study the asymptotics 
of the /c-instanton solution of {un,k) the Painleve I equation. The doubly indexed sequence {un,k) for 
n, fc = 0, 1, 2, . . . can be defined concretely by the following quadratic recursion relation: for k = 1 we have 

(1-7) = ^^^|l2^u/,iU(„+i_i)/2,o - ^(2n- l)^u„_i,i|, wo.i = 1. 

1=0 

while for general fc > 2 we have 

^ n—3 k—1 n 

""■'^ ^ To7i:2 r\ I X! Ul.kU{n^l)/2.0 + 6 Ui^rnUn-l.k-m 

(1.8) ^^^'^ ^ ^' ^ 1=0 m=ll=0 

25 25 "I 

- ^(2n + fc - 4)^w„_2,fc - 15 "^^(^ + 2n- 3)u„_i,fc|. 

It is understood here that u„/2,o = if n is not an even integer. 

The above recursion defines m„ ^ in terms of previous for n' < n or n' = n and fc' < k. The paper 

is concerned with the asymptotics of the sequence (wn.fe) for fixed fc and large n. As we shall see, resurgence 
analysis predicts that the asymptotics of (un,/c) to all orders in 1/n is given in terms of three known sequences 
{ui,k±i), i^J■l,k±l), {4>i,k-i) and two Stokes constants 5*1 and S^i] sec Equation (1.15). The first constant Si 
is known from the asymptotics of the 0-instanton (un,o) a-nd the second one appears in the asymptotics of 
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{un^k) for k > 2 and its exact value is unknown at present. Two new features appear in the asymptotics 
of {un,k) for fc > 2: the constant S-i and the presence of logn terms. These features are absent in the 
asymptotics of (w„,o) and 

Equation (1.8) defines but does not motivate the A:-instanton solution (w„.fe) to the Painlevc I equation. 
The motivation comes from the so-called trans-series solution of the Painlevc I equation. Trans-series were 
introduced and studied by Ecalle in the eighties; [EcaSla, EcaSlb, Eca85]. The trans-series solution u{z,C) 
of the Painleve I equation is a formal power series in two variables z and C that is defined as follows. 
Substitute the following expression 

(1.9) ?/(z,C) =^C^Ufc(z) 

fc>0 

into (1.1) and collect the coefficients of C'^. It follows that Uk{z) satisfy the following hierarchy of differential 
equations, which are non-linear for A; = 0, linear homogeneous for A: = 1 and linear inhomogeneous for fc > 2: 

(1-10) A 

-qu'I + Uk'Uk-k' = 0, fc > 1. 

k'=0 

Uk{z) is known as the fc-instanton solution of (1.1), and it has the following structure 

(1.11) Uk{z) = z^/^c-^^^'""^k{z), 
where A is given in (1.5) and 

(1.12) 0fe(z) = z-5fc/8^u„_,^"W4, 
Since C is arbitrary we can normalize 

(1.13) M0,1 = 1. 

This motivates our definition of (u„,fe). The trans-series (1.9) is what is called a proper trans-series, in 
the sense that all the exponentials appearing in it are small when z — > oo along the direction Arg(2;) = 0. 
Therefore, the instanton solutions Uk{z) give exponentially small corrections to the asymptotic expansion 
(1.2) and they can be used to construct actual solutions of the Painlevc I equation in certain sectors; see 
[Cos98]. 

The instanton solutions of Painleve I also have an important physical interpretation in the context of two- 
dimensional quantum gravity and non-critical string theory. As it is well-known (see for example [DFGZJ95] 
and references therein), the Painleve I equation appears in the so-called double-scaling limit of random 
matrices with a polynomial potential, and it is interpreted as the equation governing the specific heat of 
a non-critical string theory. The asymptotic expansion of the 0-instanton solution (1.2) is nothing but the 
genus or perturbative expansion of this string theory, and z~^/^ is interpreted as the string coupling constant. 
The instanton solutions to Painleve I correspond to non-perturbativc corrections to this expansion. In the 
context of matrix models, they can be interpreted as the double-scaling limit of matrix model instantons, 
which are obtained by eigenvalue tunneling [Dav93, MSW08]. In the context of non-critical string theory, 
they are due to a special type of D-branes called ZZ branes [Mar, AKK03]. The asymptotic behavior of the 
fc-instanton solution is then important in order to understand the full non-perturbative structure of these 
theories. 

1.3. The predictions of resurgence for the fc-instanton asymptotics. Our paper consists of three 
parts. 

(a) A proof of the all-orders asymptotics of the and 1-instantons of Painlevc I, using the Riemann- 
Hilbert approach. 

(b) A resurgence analysis of the fc-instantons of Painlevc I and their asymptotics to all orders. 

(c) A numerical confirmation of the predictions of the resurgence analysis. 
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In this section we state the predictions of resurgence analysis for the asyniptotics of for fixed k and 

large n. Recall from Equation (1.4) that the 1/n'' asymptotics of (un,o) involve the 1-instanton m^i for 
£ < k, and a Stokes constant ^i. Likewise, the asyniptotics of (wn.fe) for fixed k > 1 and large n involves two 
auxilliary doubly-indexed sequences (nn,k) and {i'n,k) and two Stokes constants Si and 5-1. The sequence 
{i^n.k) is proportional to the instanton sequence 



(1.14) Vn,k 

The sequence {^n.k) is defined by 

= (-l)"it„,i 



^kUn,k k>0 

fc = 0. 



n—2 2n 



(=0 1=0 

M2n+1,2 = 

ri-2 



192'"" '' '^2(n-l),2 



8 r 25 



25^(n+l)[ 64 



r 25 

|-^(271 + l)^/i„-l,3 + 12 ^ W(n+1-0/2,OM;,3 



2 ri+1 



for fc = 1, 2, 3, and by 



25, , 25A 

+ 12 Mn+l-i,m/ii.3-m + J^V''^ + l)*^",! H 



n — 3 A;— 1 n 

|l2 ^ fJ-LkU(n-i) /2.0 + 12 X! X! W,™'""-'.^-™ ^ 64^^"" ^ ^ ^ 4)^/^„_2,fc 



12(fc- l)(fc-3) - 

^ ' i=0 m=l 1=0 

25 25 25A ~i 

-— ^fc(fc + 2n - 3)^j„_i,fe + — (fc + 2?i - 4)i/„_i fc_2 + — ^(fc - 2)i^„.fc_2 \ 
16 16 8 J 

for fc > 4. We apologize for the lengthy formulas that define the doubly-indexed sequences {fJ.n,k) and 

{vn,k)- As in Section 1.2, there is a simple resurgence analysis explanation of these sequences given in detail 

in Section 4. An analysis based on trans-series solutions and resurgent properties suggests the following 

asymptotic result for the coefficients Un k, for fixed fc > 1, in the limit n — > oo: 

(1.15) 

Un,k ~« A "+^/^ — r(?l - 1/2) <^ fc + l)uo,fc+l + -1 "A*0,A,-H1 + V =-T-^ I 



+ (-i)"(fc - i)A-"-/3|-ir(n + 1/2) |.o,.-i + E n'Tlwf J 

"^^1 I ,_i n^-i n-h 1/2- m) J 



It n:„=i(n + l/2-m). 
^1 T./^ , 1 /o^ , i^i,fe-l(-^) 



9 r 

(-1)" _Lr(n + 1/2) (log n - log A) ^.o,fc-i + E 



rn™=i("+i/2 



- (-l)M— V.|Lr(. + 1/2) |(^(n + 1/2) - logn).o,.-i + £ ' 

1=1 nm=i(" + i/2-™) J 

where 

^'■^^m 

is the logarithmic derivative of the T function; see [01v97]. In the above Equation, the convention is that 
Un,k = Vn,k = for fc = 0. 

The above formula gives an asymptotic expansion for Un^k for n large, involving the coefficients 
fii^k±i and i'i,k±i. Up to the overall factor A^"r{n + 1/2), this expansion involves terms of the form l/n' 
and logn/n'. In order to get the m-th first terms of this asymptotic expansion for Un,k, we only need to 
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know the cocfRcicnts ui k±i, /i;,fc±i, i^i^k±i with I up to to. In particular, (1.15) gives an efficient method to 
obtain the asymptotic behavior of the Un,k for fixed k at large n. 

As a concrete and important example which also clarifies the above remarks, let us look at fc = 1. In this 
case, only the first line of Equation (1.15) contributes and resurgence analysis predicts that for even n, we 
have 



(1.16) U2n,i-A-'-+'/'p:T{2n-l/2){ 



. ^ , >^ (2zx,,, + (-l)V,2)^' , 



2^^^ I trn!n=i(2"-l/2-m) 

while for odd n wc have 



oo 



OO. 



(1.17) .2„+M - ^-2- V2|Lr(2n + 1/2) (2«o,2 - Mo,2 + E if"''' ' ^''^ ^"'^^^ 

This prediction will be proved in section 3 by using the Riemann-Hilbert approach. The coefficients mo,2 
and /io,2 can be easily calculated from the recursion (1.8): 

(1-18) "0,2 = ^, Mo, 2 = -1- 

6 

In fact one can obtain closed formulae for wo.fe and fj.o,k for all k (sec (4.34) and (4.36), respectively). The 
explicit result (1.18) gives the leading asymptotics, 

(1.19) - - (-1)") A-"+i/^|j.r(n- 1/2), 

Concrete examples of the implications of (1.15) for the asymptotics of the u^ k are given in section 5, where 
the formula is tested numerically. 

For some partial results on the Painlcvc I equation and the asymptotics of the coefficients of the 0- 
instanton solution, sec also [Cos98, CCOl, CK99]. To the best of our knowledge, a rigorous computation 
of the Stokes constant 5*1 has only been achieved via the Riemann-Hilbert approach [Kap04] or its earlier 
version - the isomonodromy method [Kap88, TakOO]. If at all possible, a computation of the constants Si 
and 5-1 using resurgence analysis would be very interesting. 

1.4. The asymptotics of the 1-instanton solution via the Riemann-Hilbert method. Recall from 
Equations (1.10) and (1.11) that the generating series ui{z) of the 1-instanton (wn,i) solution to Painlcvc I 
satisfies the differential equation 

(1.20) u'i'(z) = 12ui(z)uo(z) 

where uo{z) is the generating series of the 0-instanton solution (m„,o) solution to Painleve I. There are five 
explicit tritronquees meromorphic solutions of Painleve I equation asymptotic to 1x0(2;) in appropriate sectors 
in the complex plane, discussed in Section 2.1. To study the asymptotics of the 1-instanton (wn.i), consider 
the linear homogenous differential equation 

(1.21) v" = '^uv 

where 7 is a constant and w is a tritronquce solution to Painleve I. It is easy to see that Equation (1.21) has 

4 



two linearly independent formal power series solutions vf of the form 



(1.22) vfiz) = ^-i/8e±4^-''' g ^t(7)^-'"/' 

n=0 

where = ^^(7) is given by 



(1-23) = ^{t^(« - - f^) T ^ E u™,o6„^,_2™ I bt = 1. 

^ ^ m— 1 

Note that 6^(7) is a polynomial in 7^/^ and fe~(12) = The next theorem gives the asymptotic expansion 
of b~ when 7 > 3. Let 
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(1.24) B = ^V7- 

Theorem 1.1. For n large and 7 > 3, w;e have: 

The proof of the theorem uses expUcitly the 5 tritronquees sohitions of Painleve I equation and the 
asymptotic expansion of their difference in 5 sectors of the complex plane, computed by the Riemann- 
Hilbert approach to Painleve I. An analytic novelty of Theorem 1.1 is the rigorous computation of a Stokes 
constant which is independent of 7, when 7 > 3. 

Remark 1.1. There are two special values of 7 in Equation (1.21). When 7 = 12 and w is a solution of 
Equation (1.21), it follows that v is the 1-instanton solution of Painleve I and Theorem 1.1 implies Equations 
(1.16) and (1-17). When 7 = 3/4 and w is a solution of (1.21), then v = —2v'/v is a solution of the Riccati 
equation 

(1.26) 2i)' -v^ + 3u = 

studied in [GMIO]. 

1.5. Acknowledgements. M.M. would like to thank Jean-Pierre Eckmann, Albrecht Klemm, Sara Pas- 
quetti, Pavel Putrov, Leonardo Rastclli, Marco Ranch, and Ricardo Schiappa for discussions. 

2. ASYMPTOTICS OF THE 0-INSTANTON OF PAINLEVE I 

In this section we review and extend the calculation of the asymptotics of the 0-instanton sequence (unfi) 
presented in [Kap04]: 

(a) Use the 5 tritronquce solutions to Painleve I and their analytic properties obtained by the Ricmann- 
Hilbert method, as an input. 

(b) No two tritronquee solutions are equal on a sector, however their difference is exponentially small. 
Using the 5 tritronquee solutions and their differences, define a piece-wise analytic function in the 
complex plane which has (1.2) as its uniform asymptotics in the neighborhood of infinity. 

(c) Apply a mock version of the Cauchy integral formula to obtain an exact integral formula for u„.o in 
terms of the jumps of the piece-wise analytic function defined above. The knowledge of the explicit 
value of the relevant Stokes' multiplier (available due to the Ricmann-Hilbcrt analysis) yields the 
large n asymptotics of u„fl. In fact, we extend the result of [Kap04] and we obtain the asymptotics 
to all orders in 1/n. 

Note that this method consists of working entirely in the z-plane (and not in the Borel plane) , in all sectors 
simultaneously, and our glued function is only piece-wise analytic. This method is different from the method 
of Borel transforms analyzed in detail in [CK99; CCOl]. 

2.1. A review of the tritronquee solutions of Painleve I. In this section, to simplify our notation, we 
will replace u„_o by a„, and use the symbol Un riot for the n-th term in the formal trans-series expansion, 
but for certain exact solutions of the first Painleve equation specified below. 

Recall that equation PI (1.1) admits a formal 0-paramcter solution with the power expansion Uf{z) (cf. 
(1-2)), 

°° r 25n2 - 1 1 " 

(2.1) Uf{z) = z^/^ ^ a„z"''"/^ ao = 1, Oi = -1/48, a„+i = — a„ ^ 2 ^ a„ian+i-m- 

n—O m— 1 

Remark 2.1. In [Kap04], the Painleve I was studied in the following form: 

(2.2) y,, = 6y2 + x. 

The change of variables u ~ &/^y and z = e^"6~^/^x transforms (1.1) to (2.2). Comparing with [Kap04], 
we also find more convenient to modify the indexes for the tritronquce solutions. 
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Using the Riemann-Hilbert approach (sec [FIKN06, Ch.5.2] and also [Kap04]), it can be shown that there 
exist five different genuine mcroniorphic solutions of (1.1) with the asymptotic power expansion (2.1) in one 
of the sectors of the z-complex plane of opening Stt/S, see Fig. 1: 



uo{z) - Uf{z), z^oo, argze(-^,^), 

(2.3) Uk{z) = e-'^^"uo(e-'^'2) ^w/(^), ^ ^ oo, arg z € (-f + f fc, f + f fc) , 

(2.4) Uk+5{z) = Uk{z), 

where in u/(z) appearing in (2.3) the branches of z^/^ and z-S"/^ are defined according to the rule, 

(2.5) zi/2 = ^e^^g^/2, argze (-f + f fc,f + f fc), z-5»/2 ^ (^1/2) . 

The asymptotic formula (2.3) is understood in the usual sense. That is, for every e > and natural number 
N there exist positive constants Cn,^ and RN,e (depending on N and e only) such that, 



N 



(2.6) 



Uk{z) 



.1/2 



-5n/2 



< Cat., 



n=0 



5 



47r jL 



Estimate (2.6) is proven in 



Vz : |z| > i?7v,e, argz e [ ^ , ^ , ^ 

Here, the branches of z^/^ and z^^"/^ are again defined according to (2.5) 
[Kap04]. 

Remark 2.2. We note that estimate (2.6) implies, in particular, that each tritronquee solution Uk{z) might 



have in the sector [- 

i?o,e- We also note that the number Rm,^ in estimate (2.6) can be replaced by i?o,e for every TV. 



^■^-A; + e, ^ + — e] only finitely many poles which lay inside the circle of radius 




Figure 1. The sectors of the z-complex plane where the tritronquee solutions (2.3) m_2, uo, ui, U2 

are represented by the formal series Uf.\n the dotted sectors, the asymptotics at infinity of the tritronquee 
solutions is elliptic. 

Furthermore, in [Kap04] it is shown that the exponential small difference between the tritronquee solutions 
Uk{z) and Uk+i{z), within the common sector where they have the identical asymptotics Uf{z) in all orders, 
admits the following explicit asymptotic description: 

(2.7) z^oo, argze(-f,f): u,{z) - uoiz) = z-'/'e-'^^''^^! + 0{z-'/*)) 

z^oo, argze (-f + ^fc,f + ^A:): 
Uk^.iz) uk{z) = e'i 01 ^ o{z-''')) . 

Here again the asymptotic relations mean that the differences Uk+i{z) — Uk{z) admit the representation, 

(2.8) u.+i(z) - uk{z) = ec+D^^ ,-i/8,(-i)-^^.^/^ + ^ 
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whith the error term r{z) satisfying the estimate (cf. (2.6), 

(2.9) \r{z)\ < a\zr'^\ 

Vz: \z\>R,, argze [-2- + i^k + e, ^ + fk - e], 

with the positive constants Ce and Rf depending on e only. In fact, we can take Re — i?o,c- 

When dealing with the tritronquce solutions Uk{x), it is convenient to assume that k can take any integer 
value, simultaneously remembering that Uk+5{x) — Uk{x), see (2.3). In other words, in the notation Uk{x) 
we shall assume that 

(2.10) kel, mod 5, 
unless k is particularly specified, as in (2.16) below. 

Remark 2.3. The existence of the tritronquce solutions Uk{x) can be proved without use of the Riemann- 
Hilbert method (see e.g. [JKOl] and references therein). Indeed, these solutions had already been known to 
Boutroux. The Riemann-Hilbert method is needed for the exact evaluation of the pre-exponent numerical 
coefficient in the jump-relations (2.7), i.e. for an explicit description of the quasi-linear Stokes' phenomenon 
exhibited by the first Painleve equation. 

2.2. Gluing the five tritronquee solutions together. The main technical part of the approach of [Kap04] 
to the 0-instanton asymptotics is a piece-wise meromorphic function with uniform asymptotics (2.1) at 
infinity. This function is constructed from the collection of the tritronquece solutions Uk{z), k = 0,±1,±2, 
as follows. 

First, observe that the change of independent variable z = <^ turns (2.1) into the non-branching series. 



(2.11) uf{t) :=u/(t2) = ^a„t-5"+i. 
Multiplying (2.11) by t^^"^, we obtain the formal series 

oo JV 

(2.12) uf\t) = uf{t)t^«-^ =P,N-i{t)+aNt-^^ ant-'^"-''^~\ P5N-l{t) = Y.aN-J 

n=N+l ra=l 

Let us perform the similar operation with the solutions Uk{z) 

(2.13) Uk{t)^Uk{t^), 

(2.14) 4^)(t) = uu{t)t''^-^ - P,N-i{t). 



Observe that. 



u 



i''\t) = aNt-^+ON{t-''), t^oo, argte (-f + f fc,f -f f fc) 



The symbol On indicates that the relevant positive constants in the estimate depend on N. More precisely, 
estimate (2.6) implies that 



(2.15) u[^\t)-aNt-^ <Cn,M 



-6 



yt : \t\ > Rl'l argt G [-f + f fc + e, | + f fc _ e] . 

Moreover, in view of Remark 2.2, we conclude that each function u'j^\t) might have in the sector [— ^ + 
^fc -I- e, ^ -I- ^fc — e] only finitely many poles whose number does not depend on iV, and they all lay inside 

1/2 

the circle of radius i?g ^ . Indeed, by (2.14), the possible poles must coincide, for every A'', with the possible 
poles of the corresponding function Uk{t). 
Put 

(2.16) u^'^\t) = uf\t), argt e (-^ + ^fc,^fc), fc = 0,±l,±2. 
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These equations determine u^'^^{t) as a sectorially meromorphic function u^^''{t) discontinuous across the 
rays argi = ^fc, k = 0, ±1, ±2, see Fig. 2. Let us choose and then fix the parameter e = eg in such a way that 
the closure of each sector (— ^ + ^fc, ^A:) is included in the corresponding sector [— ^ + ^fc+e, ^ + ^fc— e] . 
Then, the function u*-^-* (t) would have no more than finitely many poles in the whole complex plane, the 
number of poles would be the same for all N, and they all lay inside the circle of radius 

(2.17) i?o - Ro(l 

In addition, the following uniform asymptotics at infinity takes place, 

(2.18) u'^^\t) ^GNt-^ +ON{t^^), t^oo. 

The exact meaning of this estimate is the existence of the positive constant Cn{= C'jv.eo) such that, 



(2.19) u(^)(t)-aAr<-i 



<CN\t 



Vt : 1^1 > Rq. 

Here, Rq is defined in (2.17), and it does not depend on A^. We remind also the reader that the circle of 
radius Rq contains all the possible poles of u^^\t). 

The jumps of the function u^^)(<) across the rays argi = ^fc, k — 0, ±1,±2, oriented towards infinity, 
arc described by the equations, 

(2.20) uf\t)-u^^\t) = U^^\t), argi=ffc, fc = 0,±l,±2, 

where u^^\t) and u'^''(i) are the limits of ii'^^^t) as we approach the rays from the left and from the right, 
respectively. (We note that UQ^\t) = u[^\t), as it follows from (2.3)). By virtue of (2.7), the jump functions 
J7^^-'(t) satisfy the estimate, 

(2.21) ?7W(i)| 2. =e't(^+i)^t^^-9/V-i)''+'^*'''Yl + 0(t-5/2)V 

'argt= — 2y/Tr V ^ ^7 

which, again, means (cf. (2.8) and (2.9)) that 

(2.22) m''\t)\ 2. =e^f('=+i)^f5^"9/V-i)'^'^*''Yl + r(t))y 

'argt= — V '/ 

with 

\r{t)\<C\t\-^'\ Vi:|t|>i?o, argi=ffc, 
where the constant C(= C^q), similar to i?0; is a numerical constant not depending on N . 
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2.3. An integral formula for the 0-instanton coefficients and their asymptotic expansion. Con- 
sider the integral, 

of the function u^^'> along the circle of radius R centered at the origin and counter-clockwise oriented. For 
sufficiently large R we can apply to the integrand estimate (2.19). Noticing that the constant is the 
same along the whole circle, we conclude that 



(2.23) 



that is. 



(2.24) 



— i u<-^\t)dt-aN 



< CjqR 



-5 



Vi? > i?o. 



On the other hand, since m*^^^ [t) can have only a finite number of poles lying inside of the circle with the 
radius i?Oi the circular contour of integration can be deformed to the sum of the circle of smaller radius 
\t\ = p > i?o, stih containing inside aU the possible poles of 'iS^\t), and positive and negative sides of 
segments of the rays argt = ^fc, see Fig. 3. 



Figure 3. Deformation of the contour of integration for computation of the 0-instanton W-large asymptotics 

In other words, taking into account (2.20), we have that 
(2.25) aN = ^.<l u^^\t)dt + ON{R'^) = 



= ~T^i2fV U^''\t)dt+^l u^''\t)dt + ON{R-') = 

= --Ly2f[J U^''Ht)dt + ^(f (u^^Ht)+P,N-i{t))dt + Or,{R-% 

where we notice that the integral of the polynomial P5Ar_i(t) can be indeed added to the right hand side 
since it is zero. Let us now notice that from definition (2.16) of the function vS'^\t) it follows that 

(2.26) uW(0 + P^N^t) = u{t)t'''~\ 
where u{t) is defined by the equations (cf. (2.16)), 

(2.27) uit)^Ukit), argt e (-f + f fc,f fc), fc-0,±l,±2. 

By exactly the same reasons as in the case of estimate (2.19), we conclude from (2.6) that 

(2.28) |u(t)| < C(i)|t|, \ft:\t\>Ro, 
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with a numerical constant C'^-* this time^ independent of N. Estimate (2.28), together with (2.26) imphes 
the inequality, 

|ziW(i) + P5A.-l(t)||t|=p<CWp'^-\ 

where p is assumed to be a fixed positive number satisfying, 

Ro< p< R. 

Therefore, formula (2.25) can be transformed into the formula, 

(2.29) aN ^ -^J\'^'''^ it) dt + r j^^) [p) + [R) , 

where the error terms r^^ (p) and rj^' [R) satisfy the estimates, 

(2.30) |r;^^(p)| <C(i)p5^, ViV>l, 
and 

(2.31) \rf{R)\<CNR-\ yR>Ro, 

respectively. We remind the reader that the constants Cn depends on TV only while the constant C^^^ is a 
numerical constant independent of N. In derivation of (2.29) we also took into account symmetries (2.3) 
which allowed us to replace the sum of integrals from (2.25) by a single integral. 

From (2.22) it follows that the integrand in (2.29) satisfies the asymptotic equation. 



(2.32) ;7(^) it) I = e^f C^+i) 01 ,5^-9/4^- ^ ^^^^^ ^ 

with 



|r(t)| < C^"^/^ yt:t> Ro. 

In particular, this means that the integral of C/'^^(t) along the half line [p, oo) converges. The coefficient 
and the term A^'>{p) in (2.29) do not depend^ on R while the term in view of (2.31), vanishes as 
R — > oo. Therefore, sending i? — >■ oo in (2.29) (and keeping N fixed) we arrive at the equation, 

(2.33) aN = / U^''\t) dt + rW(p), 



^The constant C(i) can be taken equal to 



Indeed, we have that 



'■<'Hp) = (u^^'Ht) + PsN-iit)) dt. 
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which in turn can be transformed as fohows. 

r /"OO 

(2.34) aN = -^J U^''\t)dt + r^j^\p) = 



5 
2Td 



e'2 



31/4 



2tti J„ \ ^' 20F 



5 ■ 31/4 

' 4^3/2 



5 ■ 31/4 
47r3/2 



^5iV-9/4^ 



8^/3 j5/2 



dt 



e'2 



2V^ 



^5W-9/4g-5^t= 



The first integral, /i, in the last equation reduces to the Gamma-function integral, 



(2.35) h 



5 ■ 31/4 
" 47r3/2 



' .5JV-9/4^-»|St^/2 _ 3^/^ / 8^/3 ^ -2Af+l/2 

2^3/2 1 5 ; 



s2^-3/2e-^ds 



3I/4 ,8\/3,-2iV + l/2 



27r3/2 



r(2iV- i). 



For the second integral, I2, we have 



(2.36) 1/2 1 



5 • 31/4 



47r3/2 



5JV-9/4 -^t'^/^ 



dt 



< 



5.3I/4 rp 



4^3/2 

5 ■ 3^/4 fP 
47r3/2 



^5Ar-9/4^^^ 



31/4 



5Ar-5/4 ^(2)P; 

47r3/2(iV- 1/4)'^ iV 



For the third integral, I3, we use the asymptotics (2.32), 

r 0I/4 /-oo 

(2.37) -^'^ ' 



r q1/4 /"OO 5A^ 

' ^' 47ry^ 7^ ^ 2> ^ ' 

where the constants C^^\ j = 2, 3, 4 are positive constants which do not depend on N. Equation (2.35) and 
estimates (2.36), (2.37), and (2.30) mean that 



(2.38) 



gn = 



3I/4 8V3.-2jV + l/2 
2^3/2 1 5 ) 



r{2N - i)(l + O(iV-i)) + 0(p5^), 



as iV — 00. This, in turn, implies the following expression for the leading term of the large N asymptotics 
of the coefficient a^. 



(2.39) 



31/4 8y3.-2JV+l/2 



r{2N - ^){1 + 0{N-'^)), N- 



27r3/2^ 5 ' 

Equation (2.39) provides us with the leading term of the large n asymptotics of a„ = u„^o- In order 
to reproduce the whole series (1.4) we note that the Riemann-Hilbert analysis of the tritronquee solutions 
yields in fact the full asymptotic series expansion for the differences between the tritronquee solutions. This 
means the following extension of the estimates (2.7): 



(2.40) 



0I/4 °° 
..+i(.) - u,iz) ^ e'f (^+1)1- .-V8,(-i)-^^.^/^ 1 + ^ 



CjiZ 



-5n/4 



00, arg z £ 



2ir I 4tt i„ 27r , 47r 



k), 
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with some coefficients c„. In its turn, asymptotics (2.40) impHes that 

0I/4 ™ 

(2.41) m^Ht)\ 2. =e»f(''^+l)±— t5^-9/V-l)'+'^*''Yl + y C„i-W2+0(t-5(rn+l)/2)y 

^ n— 1 

Therefore, the integral of [/'■^^(t) along the half-line (p, 00) can be estimated ^ as follows: 
(2.42) 

= C„A-2^^+"+l/2 e-42^-"-3/2^^^0/^-2W+™+3/2 r g-t^2W-m-l-3/2^A ] 



= |^f]^c„A-2A^+"+i/2r(2iV - n - i) + o(^-2^+'"+3/2r(2iV - m - §)) j + ©(p^^-^/^) 

Here, m is arbitrary but fixed and all the constants in the error terms depend on m and p (which is also 
fixed) only. We also have used our standard notations, 

, 8V3 . 3i 

A=^—, bi = -\—=, Co = 1. 
5 Syti" 

Substituting (2.42) into (2.33) we arrive at the final estimate for ajq-, 

To complete the proof of the - instanton asymptotics (1.4) we need to identify the coefficients c„ as the 1 
-instanton coefficients Un,i- To this end, it is enough to notice that the difference, v(z) = Uk+i{z) — Uk{z), 
of the tritronquee solutions satisfies the linear differential equation, 

(2.44) v" = 6viuk+i+Uk). 

Since both and Uk+i have the same power series expansion in the relevant sectors, we conclude that the 
formal power series solution of equation (2.44) must coincide with the formal power-series solution of the 1 
-instanton equation (1.20). Hence the desired equation, 

(2.45) c„ = M„_i, Vn, 
which completes the proof of (1.4). 

3. Asymptotics of the 1-instanton of Painleve I 

3.1. A 5-tuple of differential equations. Let us consider the first instanton term Cui in the formal series 
(1.9). This term satisfies the linear homogeneous second order ODE 

(3.1) u'l = 12uoMi, 

where uq is the 0-instanton term in the trans-series expansion (1.9). This suggests the problem of studying 
the linear equation 

(3.2) v" = -fuv, 



We omit the routine technical details which are similar to the ones carefully presented in the derivation of (2.39 ) 
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where u is one of the tritronquce solutions of the first Painleve equation (1.1) and 7 is an arbitrary positive 
constant. Indeed, we have five different Schrodingcr equations with meromorphic potentials ^Uj having the 
same formal power expansion (2.1) in the relevant sectors, see (2.3): 



-5n/2 

1^ 



(3.3) v'- = jUjVj, 7"/(^) = 7^^^^ X! ""^ 

n=0 

Equation (3.2) has two linearly independent formal solutions, 

00 

(3.4) vf = e'^J:^, =B^z^'^~\\nz, E± = ^ &±z-W4^ 



To simplify our notation, in this section, we use the symbol b~ instead of u^^i- 

Since we have five equations with meromorphic potentials, we have five pairs of solutions with the as- 
ymptotic expansions (3.4) as z — > 00. Each pair has this asymptotic expansion in the sector where the 
corresponding potential has the power-like asymptotics (2.1). The sectors are depicted in Fig. 1. Denote the 
solutions corresponding to the potential Uj{z) by the symbol v^{z). We have, 

(3.5) vf{z)^vf{z), argze (-f + f j,f + f j). 
We also observe that, 

(3.6) .f(.) = e-^^^.-^-)(e-^^.), .,(±) = {"^ "^^^ 

\T, J is odd, 

and 

(3.7) vf^^^{z)^vf{z). 

Similar to convention (2.10) concerning the subscribe k in Uk{x), we shall assume that 

(3.8) j G Z, mod 20, 

in the notation u^(z), unless j is particularly specified, as in (3.19) below. 

3.2. Gluing the solutions together. We observe that within the common sectors of the same asymptotic 
behavior, the solutions v^{z) and v^_i{z) have identical asymptotics v^{z) in all orders. Thus it is natural 
to ask: what is the difference between these two solutions of different equations? Let us introduce the 
differences 

(3.9) wfiz)=vf_^,iz)~vfiz). 

Because of (3.6), these differences are related to each other by the following rotational symmetry, 

(3.10) wf{z) = e~'ro3yj-^^^^\e-'¥Jz). 
In addition, using (3.3), w^{z) satisfy non-homogeneous linear ODE, 

(3.11) iwf^i)zz = iujwf_^ + -f{uj - Uj^l)vf_^. 
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The homogeneous part of this equation is the same as in (3.3) In the relevant sector, the non-homogeneity 
contains an exponentially small factor given in (2.7): 



(3.12) z-^co, 



argze(-f + + 



, 31/4 



2v^ 



l + 0(z 



-5/4^ 



9i = \nz, 

J J Q ' 



The general solution to (3.11) is given by the integral formula, 



(3.13) wf_,{z) = 



- C+V+ 



2^7, 



(z) - vf{z)v^ {x))-/{uj{x) - Uj-i{x))vf_^{x) dx 



where c"*", c~ and zq are arbitrary constants. However, our solution (3.9) is not general as being a difference 
between two functions with the same asymptotic expansion in certain sector. Furthermore, in the sector 
argz G (— ^ + + see (3.12), for j odd, the solution Vj'{z) is dominant and vj (z) is recessive, 

while for j even, v^' (z) is recessive and v~ {z) is dominant. Thus w^_i{z) admits the following representations: 



(3.14a) j odd: v^{z) dominant, Vj (z) recessive as argz € (— ^ + ^j, — ^ + ^j) : 

Wj-i{z) = cqv^ (z) - 2^^J '^^^ / .^u-i, ("^■'^■^^ ~ ^j-i(a;))w/(a;)w+_i(a;) dx 

If 

+ 2\/T'^'/(^) / " Uj-i[x))v^ (.t)u+_^(.t) dx. 



w^_^{z) = --^ l42L^^_,^ {vj'^^'>^'j (a;))(wj(a;) - Wj-i(x))v^_i(a;)da;, 



(3.14b) j even: wj^(z) recessive, w- (z) dominant as argz G (— ^ + ^j, — ^ + ^j) 



Wj_i{z) ^-\Vl («j^(2;)«7(^)-f'j^(2:)^'7(a;))(^tj(a;)-%-i(a;))^'jli(2;)rf2;. 



5 

Wj_i{z) = cov+{z) ~ ^y/jv-{z) I ^4^^^_^^^{uj{x) ~ Wj_i(a;))w+(a;)v7_^(x)da; 



5 

1 



+ oV^^/(^) / izL,,,, K(a;) - Wj-i(x))v^. {x)v^_^{x)dx, 

Je' 5 " 'zo 



where e' 5 '■^^ ■'zq is a finite point within the indicated sector and cq is a Stokes constant which both can 
not be determined immediately. 

When z — > 00 as argz S (— ^ + ^j, — ^ + ^j), the leading contribution to the asymptotic behavior 
of w^_-^[z) is given by z-end point of the integral term and possibly by the non-integral term. By standard 
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arguments, 



(3.15a) j odd: v^{z) dominant, (z) recessive as argz e (— ^ + ^j, — ^ + ^j) 



57^yl(-A + 2B) 
2 • 31/4 



Icie 



5/4 



1(2;) = 



3^/^ 25^7 
5^7? + 2B 



■e-(^+^)-''%-3/4(l + 0(z-5/4)), 



(3.15b) j even: v^{z) recessive, Vj (z) dominant as argz g (— ^ + ^j, — ^ + ^j) : 



""^"-1^^" 5^/^A(A + 2i3)'^ (i + C^l2 jj 



^,^,2_3^ 5,5/4^,/ + 0(;^-^/«)), i? = A/2, 

[cie^-'''z-i/8(i ^ 0(z-5/4)), < S < A/2, 

where Ci is an unknown Stokes constant. 

Taking into account the presence of the factor in the asymptotics of the functions v^{^), we define 
auxihary functions vf{^) by the relations (cf. (2.13)), 

(3.16) z = e, vf{0=vf{z)e-'^^^\ 
SO that (see (3.5)) 

oo 

(3.17) «f(0~S±(e»)-^6±r'°", e^oo, arg^e (-If + |. + JL^^-). 

n=0 

Using (3.6), we observe that 

(3.18) vfiO = e-^^oH-^(^\e-^^^^), vf_,,,{0 = vfiO- 

Let us introduce the sectorially meromorphic function V^(^^(^) discontinuous across the rays arg^ = j^j, 



j = -9,-8,...,9,10: 

(3.19) V^''HO = yJiOe''''-'-ProN-iiO, arg^e (-flj + ^j,^j)' J = -9, -8,..., 9, 10, 
where polynomial ^'ioAf-i(0 of degree lOA^ — 1 is defined by 

7V-1 

(3.20) PioN-iio--^T.^ne'^''~"^~'- 



ASYMPTOTICS OF THE INSTANTONS OF PAINLEVE I 

According to (3.17), this function has the uniform asymptotics at infinity, 



17 



(3.21) 



— !• OO, 



and, due to (3.9) and (3.16), has the following jumps across the rays arg^ = ^{j — 1), j = —9, —8, . . . , 9, 10, 
oriented towards infinity: 



jlOW-l 



(3.22) arge = - 1)^ K 'iO - ^^"^0 = {iJiO - 



5,-7, ...,11. 



(We note that i'ii(^) = v_g{z), in virtue of (3.7)). The jumps for even and odd j-s are different. Namely, 
as ^ — > CO, using (3.15), we find 



(3.23a) ^->oo, argC = ^(j-l), j odd 



A{A + 2Bf 



(3.23b) C^oo, arg^ = f^r(j - 1), j even 



v:(^)(^)_yr'(0 = <! 



3^/" 2i?V7 ^^:o 



A{A-2B] 
,.2 -31/4 



e^«^i°^-«(l + 0(r'")), i3>A/2, 



^,.^.^o ^^^.u^io^^4 + 0(r^)), i? = A/2, 



2B?"'^ioJV-i + 0(r '")) , < B < A/2. 



Remark 3.1. In this section, as well as in the rest of the paper, we skip the detail explanation of the 
standard meaning of the symbols 0{...) and On {■■■)■ This meaning has been spelled out in Section 2 (see 
e.g. (2.9), (2.15), (2.22), and (2.24)). We also omit, when performing asymptotic calculations, the routine 
technical details. They are similar to the ones which were carefully presented in Section 2. 

3.3. An integral formula for the 1-instanton coefflcients and their asymptotic expansion. Due 

to (3.21), the 1-instanton coefficient bj^, whose iV-large asymptotics we are looking for, is the coefficient at 
the term 1/^ of the asymptotic expansion of F(^)(C) at C = OO. Therefore, 



(3.24) 



^TT? 7|J|=_R>1 



Similarly to the 0-instanton case, collapsing this circular path of integration to a circle of a smaller radius 
1^1 — r containing all the singularities or branch points of V^'^-'(^) (the latter can appear at the poles of 



18 



STAVROS GAROUFALIDIS, ALEXANDER ITS, ANDREI KAPAEV, AND MARCOS MARINO 



Uj{^^) only), wc find 

(3.25) h-^^^l \/W(0d^ + Oiv(i?-^°) = 

fc=-5 Jre's'' 

-J_y e-*^(2fe+l) +((g-«-j(2fe+l)^)8N BC-^lOA^^^^0(^,lOiV)^0 (^-10)^ 

TZl J J, 111 J J. 

Here, we used that 

where the constant C does not depend on N . It is not difficult to find A^-large asymptotics of the coefficients 
b~^. Similar to the - instanton case considered in Section 2, we first send i? — >■ cx) in the last equation of 
(3.25) and then, using (3.15a), analyze the large N behavior of the resulting integrals: 



( 4 • 31/4 A{1 + (-1)^) - 2B{1 - (-1)^) 



(3.26) 



257r3/2 



- ^2 



7A-^-5r(7V-i)(l + 0(7V-i)), 7>3, 



^^-<!(-l)-|^V^(2B)-- 



^r(iv + -)(i + o(iv-V2)), ^ 



N 



ci — (25 



r^r(7V)(i 

:V3. 



-0(7V-i)), 



< 7 < 3. 



We recall that B — |^/7 and A — ^ ■ 

Equation (3.26), in the case 7 = 12, produces the leading term of the asymptotics of the 1-instanton 
sequence u„^i. The proof of the whole series (1.16)-(1.17) is similar to the proof of the 0-instanton asymptotic 
series (1.4) performed in the end of Section 2. 

Remark 3.2. Note, that for 7 > 3, the unknown Stokes constant ci does not contribute to the leading term 
of the large N asymptotics of b'^] indeed, it appears in the exponentially small terms only. For 7 = 3, this 
contribution appears in the subleading term of order N~^^'^. However, for < 7 < 3, the asymptotics of 6^ 
is simply proportional to Ci. This is also confirmed by numerical computations. 

3.4. Tritronquee solutions and the induced Stokes' phenomenon. The analysis performed in the 
previous subsections gives rise to the following observations concerning the Schrodinger equation, 

(3.27) v" = juv. 

Let us take two different tritronquee Painleve I functions, say the functions uo{z) and wi(z), as the po- 
tentials u{z) in equation (3.27). The corresponding solutions, Vq{z) and v^{z), though solutions to different 



ASYMPTOTICS OF THE INSTANTONS OF PAINLEVE I 



19 



equations, would have the same asymptotics in ah orders in the common sector, — 27r/5 < argz < 27r/5. 
The same is true for the pair Vq{z) and v^{z). Hence the possibihty arises of evaluating the exponentially 
small differences, v^{z) — v^{z) and (z) — (z). In other words, the existence of the tritronquee Painleve 
transcendents allows to introduce the notion of the Stokes' phenomenon for the collection of the solutions 
to the family of linear equations (3.27) consisting of five equations generated by five tritronquee functions 
as potentials u{z). We shall call this phenomenon the induced Stokes' phenomenon or, more lengthy, the 
induced Stokes' phenomenon generated by the first Painleve quasi-linear Stokes' phenomenon. It is worth 
noticing that no single equation (3.27) with u{z) = Uk{z) generates a meaningful Stokes phenomenon in the 
set of its solutions. Indeed, the structure of the solutions to equation (3.27) in the complementary sector, 
where the potential has infinitely many poles (see Figure 1), is extremely difficult to describe. We need five 
different potentials in (3.27) to cover the neighborhood of infinity by the sectors with the regular behavior 
of both potentials and the solutions. 

There is a threshold value of the parameter 7 in (3.27), namely, 7 = 3. When 7 > 3, the induced Stokes' 
phenomenon generated by the first Painleve quasi-linear Stokes phenomenon is completely described by the 
latter. Indeed, the origin of the coefficient 3^/^/2y^ in formulae (3.15a) and (3.15b) is the Stokes' constant 

which is featuring in the quasi-linear Stokes' relations (2.7) for the tritronquee solutions Uk{z). When 7 < 3, 
the quasi-linear Stokes' phenomenon for potentials of equation (3.27) does not control the induced Stokes' 
phenomenon for its solutions. The intrinsic Stokes' constant - the constant ci, begins to play a dominant 
role. 

There are at least two values of the parameter 7 which arc less than 3 but for which we believe an explicit 
description of the induced Stokes' phenomenon is possible. These values are*, 

(3.28) 7 = 2, and 7=^ 

In the case of 7 = 2, the two linear independent solutions of equation (3.27) admit the following representa- 
tion, 

(3.29) i^+(z) = *2i(0,-z6i/^), z;-(z) = 4-22(0, -z6i/5), 

where 4* (A, x) is the 2x2 matrix solution of the Riemann-Hilbert problem associated with the first Painleve 
equation. This Riemann-Hilbert problem is used in [Kap04] for evaluation of the Stokes constant 6*1 . Equa- 
tions (3.29) allow to use the same Riemann-Hilbert problem to evaluate the Stokes parameter associated 
with the functions v^{z). 

The Stokes constant corresponding to the case 7 = 3/4 was conjectured in [GMIO] on the basis of the 
relation of equation (3.27) with 7 = 3/4 to the symmetric quartic matrix model studied by Brczin-Ncuberger 
[BN91] and Harris-Martinec [HM90]. This relation suggests the existence of an alternative Riemann-Hilbert 
representation of the first Painleve transcendent which, through the formulae similar to (3.29) would generate 
solutions of equation (3.27) with 7 = 3/4. We expect this alternative Riemann-Hilbert problem can be 
deduced with the help of the relevant double scaling limit of the Lax pair for the skew-orthogonal polynomials 
associated with the symmetric quartic matrix model [BN91]. 

4. Trans-series 

4.1. Alien derivatives and resurgence equations. An alternative method to obtain the asymptotic be- 
havior of the instanton solutions to the Painleve I equation is based on the ideas of resurgence introduced 
by Ecalle. In this approach, one first constructs the trans-series solution to the differential equation and 
computes its alien derivatives. Using these, one can easily deduce the asymptotics by using contour defor- 
mation arguments. This approach has been used in [CK96, CK99, CGll], see also [GMIO], section 4, for an 

^ Painleve I equation (1.1) provides the linear equation (3.2) with the potential u{z). We note that the value of the parameter 
7 in equation (3.2) can not be changed via the scaling of the variables in the Painleve equation without violation of the chosen 
form of the latter. The values (3.28) corresponds to form (1.1) of the equation PI chosen in this paper. 
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application to a first-order differential equation of the Riccati type. In the case of Painleve I we deal with 
a resonant equation, i.e. there is an integer linear combination of its eigenvalues which is null. In this case, 
the method of resurgence is slightly more subtle. In particular, the relevant trans-series solution involves 
logarithms, as we will see. 

As we mentioned in the Introduction, the instanton solutions to Painleve I appear as trans-series solutions, 
see (1.9). In order to understand their asymptotics we need however a more general trans-series solution, 
which is a formal, fwo-parameter series of the form 

(4.1) u{z,Ci,C2)^ c'lC^y^nUz). 

n,7n>0 

The Un\„i{z) have the structure 

(4.2) u„|,„(z) = zi/2e-("— ^^(^)^ 



where 4>n\m{z) do not contain exponentials c^^^ ' . When C2 = 0, Ci = C we recover the trans-series 
solution (1.9), therefore 

(4.3) u„|o(z) w„(z). 

Notice that the more general trans-series (4.1) is not proper (in the sense explained in the Introduction), 
since for any given direction in the complex plane, there is an infinite number of terms in (4.1) involving 
exponentials which grow big as 2 ^ 00 along that direction. Therefore, the trans-series (4.1) is more general 
than those considered in for example [Cos98]. 

If we substitute (4.1) into (1.1) and collect the coefficients of C"C™, we find that the u„|m(2:) obey a set 
of coupled, inhomogeneous linear ODEs 

(4.4) - 6 ^ Wfc|/Mn-fe|m-i = 

fe,/>0 

generalizing (1.10). In the following it will be convenient to introduce the variable 

(4.5) X = z^/^ 

We now compute the alien derivatives of these solutions. The alien derivative A^^ of a formal power 
series 4>{x) was introduced by Ecallc [EcaSla, EcaSlb, Eca85]. To obtain A^^c/), one essentially computes the 
discontinuity of the Borel transform of (/)(x), 0(C), at the cut in the Borel plane starting at C = wA. This leads 
under suitable assumptions to a series in C{C}, whose inverse Borel transform is the alien derivative Ai^0, 
see also [CNP93, SS03] for more precise definitions. A crucial property is that the pointed alien derivative 

(4.6) = e"^-A^ 

commutes with the standard derivative. This makes possible to relate alien derivatives to trans-series solu- 
tions. In our case, if we apply the pointed alien derivative to the Painleve I equation we obtain the linear, 
second order ODE 

1 d2 . 

(4.7) - - — (A„u(z,Ci,C2)) + 2w(z,Ci,C2)A^m(z,Ci,C2) =0. 

o dz^ 

This equation has two linearly independent solutions: 

9k(z,Ci,C2) du{z,Gx,Gi) 

and we conclude that 

(4.9) A..(z,C„C2) = a.(C„C2) ^"^^;,^^'^^^ 6.(C„C2) ^"^^;,^^'^^^ 

where aa;(C'i,C2), &cj(C'i,C2) are in principle formal power series in C\^Ci but they are independent of z. 
This type of equation, relating the pointed alien derivative to trans-scries solutions, is called in Ecalle's 
theory the bridge equation; see [CNP93, SS03] and particularly [GSOl] for an example of a second-order 
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difference equation. In order to understand the asymptotics, we are are particularly interested in the cases 
a; = ±1. Let us first analyze the case lj ~ 1. Wc find 



J2 e-("+i-™)^^Ci"CrAi(/)„|„=ai(Ci,C2) ^ nC^'C^c-^^-'^''^^^,^ 

n,m>0 



(4.10) 

+ 6i(Ci,C2) J2 "^CrC2"'ie-("-™)^^(/)„|^. 

m>l,n>0 

Let us first look at the term multiplying c~^^ in both sides. In the left hand side this corresponds to n = m, 
and wc obtain 

^(CiC2)"Ai(/)„|„ ^ ai(Ci,C2) J2{n + l)(CiC2)"</.„+i|„ 

(4.11) 

m>l 

Since the left hand side is only a function of C1C2, we conclude that 

(4.12) ai(Ci,C2) -;^al^fc(ClC2)^ &i(Ci, C2) = C| ^ 6i,fc(CiC2)'=. 

fc>0 fe>0 

If wc now equate the different powers of c"^^ and Ci, C2, wc find the equation 

min(n,m) min(n,m — 1) 

(4.13) Ai(/)„|„ = ^ ai,fc(n + 1 - A:)0„+i_fc|m_fe + ^ fci,fc(?7i - 1 - fc)(^„_fc|„,_i_fc. 

A;=0 fc=0 

The case a; = — 1 is very similar. We have now that 

(4.14) a_l(Cl,C2) = C2^a_l,,(ClC2)^ C2) - ^ fei,;.(CiC2)'= 

*;>0 A;>0 

and 

min(n — l,m) min(n,m+l) 

(4.15) A_i</)„|^ = ^ a_i,fe(n - 1 - fc)(/i„_i_fe|™_fc + ^ 6_i,fe(m + 1 - /c)(?!)„_fc|™+i_fc. 

fc=0 k=0 

We are particularly interested in the solutions with m = 0, corresponding to the instanton solutions of 
Painleve I. In this case, the equations for the alien derivatives read (we denote 0fc|o = 4'k) 

(4.16) Ai0fc = 5i(fc + l)(/.fe+i, k>0 
where = a±i,Oj fmd 

(4.17) A_i0fe = 5_i(fc-l)0fe_i + 5_i</.fc|i, fc>0, 

where S^i = 6-1,0 a-nd we understand that = in the case k ^ 0. The three unknown constants S±i, 
S-i are the Stokes constants for the Painleve I equation. 

As we will see, a consequence of these equations is that the asymptotics of the instanton solution Uk to 
Painleve I is determined by the solution u'^., with k' = fc ± 1, and by the trans-series solutions Uk\i, which 
do not belong to the instanton sequence. 

4.2. A study of the trans-series. Wc will now study the trans-series solutions It turns out 

that there are three differerent cases: n = 0, n = 1 and n > 2. We will now study them in detail. 
For n = 0, Uo|i satisfies the same linear ODE than wi, 

(4.18) -i<|i + 2uowo|i = 0, 
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which indeed has two hnearly independent solutions: one of them, corresponding to ui, is exponentially 
decreasing along the direction argz = 0, \z\ oo. The solution corresponding to Uo|i is exponentially 
increasing along this direction, and it is given by 

(4.19) ^.0|lW-2-^/'c+^^^'Vl(^), 
where Hi{z) is a formal power series in z^^/^; 

(4.20) A^i(z) = ^/i„a^"'"/' 

Tl>0 

and we normalize again 

(4.21) fio,i = 1. 

The coefficients fi„,i satisfy the same recursion than (1.7), with the only difference that we have a 

minus sign in the right hand side. One immediately finds that 

(4.22) //„,i = (-l)"u„4. 

Let us now consider which satisfies the linear inhomogeneous ODE 

(4.23) - + 2iioWi|i + 2uiUo\i = 0. 
It is easy to see that wi|i(z) has the following structure 

(4.24) Ul|l(z) = z-3/4^;,„,2^-^"/^ 

n>0 

and in addition /i2n+i.2 = 0. The even coefficients satisfy the recursion 

n— 2 2n 2^ 

(4.25) /i2„,2 = - ^ m,2Un-lfl - ^(-l)'lt/,l"2n-;.l + J^i'^'^ ~ l)^M2(n-l),2 

and we find, for the very first terms, 

V4/ 75 300713 , 

For fc > 2, the formal solutions Ufc|i develop a new feature: they contain logarithms. This is due to the 
resonant character of the Painleve I equation. These solutions have the following form: 

5 

(4.27) Ufc|i(z) -logz.gfc_i(z) + /fe+i(z), fc > 2, 

where 

ri>0 

4.28 

n>0 

The factor 5/4 in (4.27) is introduced for convenience, since in the resurgent analysis it will be convenient to 
use the variable x in (4.5). The functions fk,gk appearing in (4.28) satisfy the coupled system of equations 



6 

(4.29) 



^ k-l 

-g'l + 2uo5fc + 2 X! ""ifffe-i = 0' fc > 1. 



- \ri + 2wo/fc + 2 uJk-^ + -^9k-2 - =0, fc > 3 
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In the second equation we set fi{z) = Uq|]^(z), f2{z) — Ui^i{z). It is easy to see, from the recursion relation 
obeyed by the coefficients i'n,k-, that 

(4.30) Vkiz) = Ckukiz), k > 1 
where C is a constant given by 

(4.31) C=i§. 

The value of this constant can be fixed by looking at the equation for ^211(2). 

Finally, one can easily find recursion relations for the coefficients /i„.fc appearing in fk{z)- The cases A: = 3 
and fc > 4 are slightly different. For /c = 3, one finds 

n-2 



(4.32) 



" 2bA{n + 1) ^ 1)'m»-i,3 + 12^ U(„+i-o/2.o/^^3 

2 "+i 25 2bA 1 

m=l (=0 



i=0 

while for fc > 4 we have 



(4.33) lZ[k^l){k 6)1 ™=u=o 

25 25 25A 1 

— Afc(fc + 2?! - 3)^„_i.fc + — (fc + 2n - 4)i/„_i,fc_2 + -^-(fc - 2)i/„.fc_2 f 
16 16 8 J 

One interesting aspect of the doubly-indexed sequences {un,k), i^J■n,k) is that we can also consider their 

dependence on k for fixed n, and it turns out that one can find closed formulae for these coefficients from 

the recursion. For example, (ito.fc) is given by 

(4.34) uo,k = (12)i-'=fc, fc > 1, 
while (ui.fe) is given by 

(4.35) = -^^--=fl09fc2 - 120fc + 24), k>2. 

16v3 ^ ' 

These formulae can be also obtained from the results in section 5.2 of [CCOl]. Using their results one can 
see that, as a function of k, 12^Un,k is a polynomial in k of degree n. Finally, we give a general formula for 
Ho^k when > 4: 

(4.36) Mo,fe = 12"''+^(A:-2)(141/c- 402), fc > 4. 

Let us close this section with a problem dealing with the physical interpretation of the full trans-series (4.1). 

Problem 4.1. As we mentioned in the Introduction, the series it„|o(2) can be interpreted in terms of the 
double-scaling limit of instantons in the matrix model, and as an amplitude associated to a ZZ brane in 
Liouville gravity. What is the interpretation of the more general trans-series m„|„(z), with m > 0, in 
the context of matrix models and in the context of non-critical strings? Do they correspond to new non- 
perturbative sectors of these theories? 

5. ASYMPTOTICS 

5.1. Asymptotics of the multi-instanton solutions. The asymptotics of the coefficients Un.k can be 
obtained from a well-known application of Cauchy's theorem in the Borel plane (see [CK99, CGll, GMIO] 
for examples). To do this, we consider the formal power series 0^ appearing in the fc- instanton solution, 
which we write in terms of the variable x = z^^^ as 

(5.1) (/.fc(x) =x-'3'=^u„,fcx-" 

n>0 
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Here, (3 — 1/2. The Borcl transform of 4)k{x) is 



n+pk-1 



Therefore, 
(5.3) 



r(/3A: + n) /q 27ri ' 

We now apply the standard deformation contour argument. The Borcl transform (pk has branch cuts at 
p = ±A, and we can deform the contour to encircle these (we also pick a vanishing contribution from a circle 
at infinity). But the contribution to the integral is given precisely by the discontinuity across the cut, i.e. 
by the alien derivatives that we calculated in (4.16) and (4.17). We then obtain the asymptotic formula 

(5-4) T^..,,. , . 7^ik + l) / dp „. . (fc-1) / dp ~ — / dp- 



r(/?fc + n) " 27ri^ ^ 7a 27ri ^ ^ p/3fc+n 27ri 7-oo pl^^+'^ ' 

For k >2 this formula involves the Borcl transform of (/)fc|i(x). To compute this we have to use 

(5.5) B(\ogxx-A^-^\ogp+pf\p'^-\ 

The calculation of the integrals appearing in (5.4) is very similar to the calculation in [GMIO]. The only new 
ingredient is the logarithm appearing in the Borel transform (5.5), which leads to an integral of the form 

Notice that the term involving ip{r + (fc — l)/3) will cancel against the contribution coming from the second 
term in (5.5). For n large we can use the asymptotic behavior 

(5.7) il;(n-r + /3)^\ogn+^^—!—^+0{l/n^). 

n 

Before presenting the general result for the asymptotics derived from (5.4), let us analyze in some detail 
the case fc = 0, since this will fix one of the Stokes constants. For k = Q the only contributions to the 
asymptotics come from ui and Uoii, and we obtain 



Un/2 



(5.8) 



„ ^ ^-"+i/2|Lr(„ - 1/2) («o.i + V ""-^ ) 



+ (-^)-"^^/^|[r(n - 1/2) {,0,1 + E . 7^^Y/2 J 



Since w„/2,o = if n is not even, the right hand side of this relation must vanish if n is odd. Using (4.22) we 
find that this is the case provided that 



(5.9) (-l)^5_i = Si 

This is exactly as in the ODE studied in [CCK04]. We then find the result [JKOl] (see also [GLM]) 

nLi(2n-l/2-m). 
Notice that 5'i has aheady been evaluated in (2.39) and it has the value 



(5.10) u„,o^A-2»+i/2^r(2n-l/2) I + Et^^ 7^^"^^. " 

1=1 
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Let US now analyze k = 2. The asympotics of w„^i involves uq and Using (5.9) wc can write it as 

and it depends on the parity of n. Using the explicit results (1.18), it is easy to check that the leading 
asymptotic behavior is as stated by Joshi and Kitaev in [JKOl, Prop. 16]. The formula (5.12) gives in 
addition the all-orders expansion of m„.i as an asymptotic series in 1/n. 

We can now write a general formula for the asymptotics of Un,k when k > 2. Using again (5.9), and 
absorbing a factor (— 1)^^/^ in S-i, we find 
(5.13) 

,-n+fiSl f , / ^^n , >fi ((fc + l)ui,k+i + (-l)"+V;,fc+i)yl' ) 
Un,k ~ A ^^—r[n - (3) i{k + 1 wo,fc+i + -1 /^o,fe+i + > — ] ^ > 



I 7-1 (n + B-m)) 



r n;„=i(^+/?-H. 

- (-l)M— ^|Lr(n + /?) (logn - log A)|.o,.-r + £ ^Z'-^^^"^^' 1 

1^ 1=1 U„r=i{n + 13 - m) ) 

- (-1)"A— ^|Lr(n + |(^(n + /3) - logn).o,.-i + £ ^^^^4^^^-^-/,.-i(-^)' j 

for k > 2, and we recall that /3 = 1/2. As compared to the asymptotics for fc = 0, 1, the asymptotics for 
k >2 involves logarithmic terms. In fact, the dominant term in the asymptotics is precisely the logn term. 
In the last line we use the asymptotics (5.7) for the ip function. Finally, notice from (4.30) that one has the 
relation (1.14) for the coefficients I'n^k- 

5.2. Asymptotics of the multi-instanton solutions: numerical evidence. We will now perform nu- 
merical tests of the predicted asymptotic behaviour (5.13) for the instanton series Un^k- The standard 
technique to do that is the method of Richardson extrapolation. This method goes as follows. Let us assume 
that a sequence s„ has the asymptotics 



(5.14) s„ 



n 

k=0 



ak_ 

k 



for n large. Its N-th Richardson transformation s[^^ can be defined recursively by 



(5.15) 



(N) _ (N-l) n_, (N-l) _ <N-1U 



The effect of this transformation is to remove subleading tails in (5.14), and 
(5.16) 

The values s^"^ give numerical approximations to oq, and these approximations become better as N , n 
increase. Once a numerical approximation to ao has been obtained, the value of ai can be estimated by 
considering the sequence n(s„ — oq), and so on. 

The method of Richardson extrapolation can be applied verbatim to verify the asymptotics of Un,i written 
down in (5.12). Let us illustrate this with a nontrivial example. According to (5.12), the sequence 

r /II A2n-l/2 r) ^ 
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is of the form (5.14) and asymptotes, as n oo, the value 

(5.18) (2u2 2 + M2 2)A^ = -— w -0.5729166666666666... 

96 




Figure 4. A plot of the sequence (5.17) and its Richardson transformations. In this and subsequent plots, 
the horizontal axis represents the integer n, and the vertical axis the values of the sequence. 

In Fig. 4 we show a plot of the sequence (5.17) and its first two Richardson transformations for n up 
to 200, which converges to the expected value (5.18). Taking n = 250 and ten Richardson transformations 
gives the numerical approximation 

(5.19) s^ol, = -0.572916666666666667... 

Notice that this test already verifies that the general trans-series solutions u„|m appearing in (4.1) arc the 
relevant objects to understand the asymptotics, since fc = 1 involves the trans-series 

Let us now study the asymptotic behavior of the instanton sequences Un^k with k > 2. The main novelty 
here is the presence of logarithms in the asymptotics, and this does not fit a priori into the standard 
framework of Richardson transformations. However, one can transform the sequence and put it in a form 
which is amenable to an analysis with standard Richardson transformations, as pointed out in [ZJ81]. Let 
us assume that we have a sequence im with the asymptotics 

(5.20) i,n ~ logms,„ -I- i„i, Sm = ^ = 

fc>0 ' A;>0 

This type of asymptote behavior appears in instanton corrections in Quantum Mechanics. The leading 
behavior of this sequence is determined by the coefficient ao , and wc would like to find a method to extract 
it numerically. To do this we consider the sequence 

(5.21) i^ = m{i„,+i-im), 
which has the asymptotics 

(5.22) im ~ log mSm+ im, Sm = X! "~T ' ^'m ^ OQ + , 

k>l k>l ^ 

It is now easy to see that, if we apply the Richardson transformation (5.15) twice to this sequence, we remove 
both the tails in l/n^ and the tails in log n/n^. This then allows a precise determination of the leading term 
flQ. Once this has been determined, we can extract the other coefficients in (5.20) by subtracting from the 
original sequence the parts of the asymptotics which are under control. 

Let us now apply this idea to the sequence of instantons of Painleve I. The first step is to consider the 
auxiliary sequence 

A"+^Un k 

(5-23) an,k = — rik 
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whose leading asymptotics is 

'27ri 



(-l)"an,fc 



+ 0{l/n') \ logn 



(5.24) 



5*1 S 1 

+ TT^ log Ai^cfc-i + TT^ik - l)uo,fe-i 

ZTTl ZTTl 



l0gAz^l_fc_ 



ZTTl I n 



+ ((^ + l)(-l)"wo,fc+i + MO,fc+l) 

+ 0{l/n^) 

Again, it depends on the parity of n, i.e. wether n — 2m or n = 2m + 1. Let us denote 
(5.25) a^|j. = a2m.k, ai°|fc = a2m+i,fc- 

In both cases their asymptotics is of the form (5.20) and we can use the method of [ZJ81] to analyze the 
sequence numerically. 

For simplicity, we will illustrate the asymptotics by focusing on the sequence with n = 2m even, and we 
will test the four leading terms displayed in (5.24), i.e. the (leading) term in log n, and the terms in log n/n, 
constant, and 1/n. In terms of the general structure written down in (5.20), we will test the values of oq, ai 
and bo,bi. These coefficients involve all the trans-series solutions appearing in the resurgent analysis. 



0.0686 
0.0685 
0.0684 
0.0683 
0.0682 
0.0681 




0.00075 
0.00070 
0.00065 
0.00060 
0.00055 
0.00050 



Figure 5. A plot of the sequence (5.26) and its Richardson transformations, for A: = 2 (left) and fc = 5 
(right). 



(5.26) 5,l^\ = s{a2m+2,k - a2m,fc) 7^{k - 1)2(12)2-'= + 0(l/m, logm/m) 



We will first test the leading term of order logn. Following [ZJ81], we first construct the sequence 

where we have used the explicit value for I'o.fe-i derived from (1.14) and (4.34). To remove tails, we perform 
Richardson transformations in the sequence (each transform is performed twice to remove both types of tails, 
as explained above). In Fig. 5 we plot it for k = 2,5, together with the second and the fourth Richardson 
transformations, and for m = 200. The convergence towards 

(5.27) - « 0.068228352037086... 

7ri\/3 

for fc = 2, and towards 

c 

(5.28) ^ w 0.00063174400034... 

IOSttIVS 

for fc = 5, is manifest in the figures. More precisely, for m = 250 and with = 10 Richardson transforma- 
tions, we obtain the numerical approximations 



(5.29) 



0250.2^"^ = 0.068228352037087.. 
^5^5°^ = 0.00063174400031... 
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Wc next test numerically the coefficient of log s/s. We consider the sequence 

(~(e) , Si 



(5.30) V,, = m(52. + -^(fc - 1)^(12)^-'=) 

whose leading asymptoties is of the form 



(5.31) 



4^ 
SttI 



(fc - l)Mi,fc_i logm ■ 



so we can apply the above procedure. The sequence hm,k for k = 2,5, and up to m = 200, together with its 
second and fourth Richardson transformations, is displayed in Fig. 6. In both cases we have convergence to 
the predicted values 

(5.32) k = 2: -0.00426427200231..., fc = 5 : -0.000847589867... 

For m = 250 and with = 10 Richardson transformations, we obtain the numerical approximations 

UJfJ = -0.00426427200235..., 

(5-33) 

"250,5 



-0.000847589866. 





Figure 6. A plot of the sequence bm,k, obtained from (5.30), together with its Richardson transformations, 
for fc = 2 (left) and fc = 5 (right). 

We can now study the constant term of the asymptoties, which also makes possible to obtain a numerical 
determination of the additional Stokes parameter S-i. We consider the sequence 



(5.34) 



C„^,k = a2m,fc + ^{^(^ - 1)2(12)^ 



2m 



log(2m) 



where we have subtracted the leading logs. According to the predictions of resurgence, as m — > oo this 
sequence asymptotes to 



(5.35) 
where 
(5.36) 



1^(^-1)^(12)2-'= 



S'_-Y - — S-^i -\ ^ — *Si. 



V3 



The numerical analysis in the case of fc = 2 gives a numerical determination of the unknown Stokes constant 
S'_i (hence, of S-i), and we find, numerically. 



(5.37) 
and consequently 
(5.38) 



27ri 
27ri 



0.31873285573864121. 



0.3882786818052856841. 



Using this value, we can verify the asymptotic behavior (5.35) for the sequence (5.34) for higher values of fc. 
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Finally, wc consider the coefficient of f/m. To do this, we construct the sequence 



(5.39) 

which should asymptote 
(5.40) 



'rn+l,k 



Si 
27ri 



T(^{k + l)uo,k+i + Mo,fe+i) - '2^^^ ~ 1) 



l)ui,/c-i 



This involves the tran-series coefficients /^o.fe+ii smd we use the numerical determination of S'_i obtained 
above. In Fig. 7 we show the sequence (5.39) for fc = 2 and A: 5, up to to = 200, together with its second 
and fourth Richardson transformations. They clearly asymptote the predictions of resurgence, 

(5.41) k = 2: -0.002295145874084..., fc = 5 : -0.0033633587118... 

For m ~ 250 and with = 10 Richardson transformations, we obtain the numerical approximations 



-0.00215 t 

-0.00220 

-0.00225 

-0.00230 

-0.00235 

-0.00240 

-0.00245 



50 100 150 200 




Figure 7. A plot of the sequence (5.39) and its Richardson transformations, for A: = 2 (left) and fc = 5 
(right). 



(5.42) 



450,2 = -0.002295145874083. 



j(10) 
^250,5 



-0.0033633587119. 



We believe these numerical tests confirm in a very clear way the predictions from the resurgent analysis 
for fc > 2. 
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